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We construct two distinct classes of exact type III solutions of the D — 4 Einstein- Yang-Mills 
system. The solutions are embeddings of the non-Abelian plane waves in spacetimes in Kundt's 
class. Reduction of the solutions to type N leads to generalized pp and Kundt waves. The geodesic 
equations are briefly discussed. 
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I. INTRODUCTION 

Exact solutions of Einstein equations have always 
attracted much attention. It is somewhat surprising 
to find exact solutions of such non- linear equations. 
Many of them were collected in the by now classic book 
by Kramer et al, which has recently been revised EI- 
Among others one finds the non-diverging spacetimes 
which were studied for the first time by Kundt in 0. 
Solutions which admit a non-expanding and twistfrcc 
null congruence are therefore said to belong to the 
Kundt's class. Such metrics are algebraically special, 
i.e., of Petrov type III, N, O, II or D. If one restricts the 
solutions to those with plane wave surfaces they are con- 
strained to be (at most) of type III; they can degenerate 
to type N and O. Conformally fiat solutions include 
for example the Bertotti-Robinson JjJ 0] universe and 
the Edgar-Ludwig metric An Einstein- Yang- Mills 
solution of type N was given long time ago in ||j and 
recovered for D — 4 supergravity in [7j. Type III solu- 
tions have somehow attracted less attention. Yet, some 
references considering the type III Einstein-Maxwell 
system can be found d, @. To our knowledge, the 
solutions presented in this note are the first type III 
solutions of Einstein gravity coupled to a Yang-Mills 
field. 

We also want to point out that the Kundt's class can be 
generalized to include a non-zero cosmological constant. 
Such solutions of typ e N were first considered in , fol- 
lowed by [H|ll El El- Type N Einstein-Maxwell solu- 
tions were generalized to an arbitrary number of dimen- 
sions in [l5| and further generalized to Lovelock- Yang- 
Mills solutions in Not very long ago the A ^ 
generalization of type III spacetimes in EI was presented 
in [13. 



II. YANG-MILLS IN KUNDT SPACETIMES 

We consider Kundt's class of metrics, described by the 
line element: 



Here, P and H are real functions while W is complex; 
u = (t — z)/y2, v = (t + z)/y2 are light-cone coordi- 
nates and z = (x + iy)/y/2, z — (x — iy)/\/2 complex 
conjugate coordinates in the transverse plane. We 
consider solutions of type III. In this case we can take 
P = I without loss of generality. As a consequence the 
Gaussian curvature of the wave surfaces u=constant, 
K = 2P 2 (ln P), Z z, vanishes. Type III solutions are 
therefore characterized by plane wave surfaces. Details 
on specific aspects of the metrics Q can be found in EI . 

It can be seen that two distinct cases arise for vacuum 
solutions of type III (for references, see El)- In the first 
case the function W does not depend on the light-cone 
coordinate v. In the second case the dependence is of 
the form W, v = — 2/(z + z). These two cases correspond 
to the Newman- Penrose quantity r being zero or ^ 0, 
respectively. 



Case W, v = 



Consider the metric QJ with: 



P 
W 



= 1 



W(u,z) 



H = 



1 /„. 



H u {u,z,z) + -(W,- Z +W, z )v 



(2a) 
(2b) 

(2c) 



Here, W is an arbitrary complex function and H is a real 
function. This is a vacuum solution of Einstein equations 
if the function H° satisfies: 



ii? 33 = H,% -Re (Wj +WW rz - z +W rzu ' 







(3) 



Numerical indices are tetrad indices 1 . All other compo- 
nents of the Ricci tensor vanish. This class of solutions 
degenerates to type N when ^ 3 = ^W, zz — 0. In this 
case one can use the remaining coordinate transforma- 
tion (see Q) to make W = 0. The metric reduces then 



ds 2 = — 2du (dv + W dz + W dz + H du) 
+ 2P~ 2 dzdz, P, v = 



(1) 1 Tetrad vectors as in 0: ei 
W dz — (H + WW) d v , e 4 = 



e 2 



du + W d z + 



2 



to that of familiar pp waves [l8| , 

ds 2 = 2 dz dz — 2 du dv — -ff (u, z, z) du 2 (4) 

where H satisfies H, zz = 0. 

The function W can have an arbitrary u dependence 
but must be at least quadratic in z so that $3 ^ 0. This 
holds for pure radiation solutions as well: they only 
enter the solution through the function H° so W is as in 
vacuum. 

We can generalize the solution to include a null Yang- 
Mills field ($ = $i = 0) of the form: 

A = a a {u 1 z, z)T a du, a a = x a (u, z) + x a (w, z) (5) 

Here, x° are arbitrary complex functions. This field was 
considered in curved spacetime for the first time in Q. 
The Yang-Mills equation in flat space is just: 



Ot, Z z = 



(6) 



This equation remains unchanged in the spacetime QJ- 
(|2"cf because of two reasons. First, the absence of com- 
ponents of the YM field strength of the form F ufJ - in this 
geometry. In the second place, the only non-vanishing 
Christoffel symbols of the form T pv p are precisely T pu p . 
As a result only the ordinary derivative survives in the 
(curved) Yang-Mills equation: 

d x F x » + T pu PF up - [A u , F UfJ/ ] = 2a, z - z = (7) 

Of course Eq. (|3J) has to be modified to account for 
the energy-momentum tensor of the Yang-Mills field: 



(8) 



Here, ^ a b is the invariant metric of the Lie group. In 
tetrad indices we have T33 = T uu . Therefore, Eq. © 
transforms to: 

H,% -Re (W 2 +WW, ZZ +W, ZU ) = 2 lab a* a* 

It is straightforward to solve for H°: 

H°(u, z, z) = f(u, z) + f(u, z) + 2 labX a X h 



Re{{W m +WW, s )z} 



(9) 



Here / is an arbitrary complex function. Eqs. JSJ and 
(|TJ)- (l2*c|) with an arbitrary (at least quadratic) W and 
the given H° describe an exact type III solution of the 
EYM system. It is a generalized Goldberg-Kerr solution. 
The type N reduction of this solution has been known 
for a long time [(J. 



B. Case W, v = -2/(z + z) 

We consider now the metric with: 
P = 1 

2v 



W = W°(u,z) 



z + z 



W° + W 
H = H°(u, z, z) + _ v - 



(10a) 
(10b) 

(10c) 



Here, W° is an arbitrary complex function and H° is a 
real function. This is a vacuum solution of Einstein equa- 
tions if the function H satisfies the differential equation: 



2^33 



(z + z) 



H° + WW 



W,° z W,°,= (11) 



This is the only non-vanishing component of the Ricci 
tensor 2 . This class of solutions degenerate to type N 
as well when $3 = W,® /(z + z) =0. In this case one 
can transform W° to zero Q and the resulting metric is 
known as Kundt waves. 

We consider the field JSJ in this background. Features 
in this geometry are the same ones as in the previous case, 
except for an extra term in the Yang-Mills equation: 



{d z + T u ™ + T vz »)F™ + (z~z) = 



(12) 



After further inspection this extra term is seen to vanish 
as T uz u = -T vz \ T us u = In this way, the 

Yang-Mills equation in the spacetime (J), ijlObfl . ll(Jc|l is 
still the same as in flat space, Eq. JBJ). Finally, Eq. (|11|) 
has to be modified to include the YM energy-momentum 
tensor T 33 = (7 q6 /27t) a,% a,|. 

We refine our ansatz in order to solve explicitly for H°. 
We choose the specific function W° — g(u)z, g(u) being 
an arbitrary complex function. This is the simplest form 
of W° such that * 3 is non- vanishing. On the other hand 
we limit the Yang-Mills field (JSJ to: 



X a = A» z 



(13) 



Here X a (u) are bounded complex functions. This field 
was considered in Minkowski space for the first time 
in [l9| and is referred to as non-Abelian plane waves 3 . 
The field l|13|) is the only form of (JSJ which gives rise to 
a bounded energy-momentum tensor. Under these as- 
sumptions the equation to solve reads: 



(z + z) 



ggzz 



g(u)g(u) = 2 lab X a (u)\ b (u) 



2 The superscript in jfaa denotes the u-independent part of -R33. 

3 The original field in 111! was \ a = \ a (u) z + some function of u 
but the later term can be gauged away without affecting gauge- 
invariant quantities, as pointed out by Coleman. 
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It is not difficult to see that: 



for x(u), v(u) can be obtained. 



H° = (f(u, z) + f(u, z)) (z + z) - ggzz 
+ a(u)(z + z) 2 {ln(z + z) - 1) 



(14) 



Here / is again an arbitrary complex function and a is 
the real function a(u) = 2j a bX a (u)X h (u) + g(u)g(u). Eqs. 
©, CSJ and O, jnSJ-lQIEJl with the given W°, H° de- 
scribe another exact type III solution of the EYM system. 
The corresponding vacuum solution (A a = 0) has to our 
knowledge not been considered before. The case g(u) = 
yields the type N reduction and can be regarded as gener- 
alized Kundt waves. A related type N solution was very 
recently given in |16| , for a Yang-Mills field of arbitrary 
higher (polynomial) dependence on z, z\ our (type N) so- 
lution is obtained by translating the line element into the 
Kundt canonical form and considering N — 0, D = 4. 



C. Geodesies 

For completeness, we present the geodesic equations 
corresponding to the two classes of solutions found, 
without attempting to provide detailed solutions here. 

• Case W, v = 

In order to get a feeling for the behavior of geodesies we 
choose the specific W = z 2 . This is the simplest W such 
that the solution presented in subsection A is of type 
III. The geodesies are described in terms of real spatial 
coordinates x, y by: 

u - V2xu 2 = (15) 

x + u 2 (H lX -x{x 2 -y 2 )) +2V2yiiy = (16) 

y + u 2 (H, y -2yx 2 ) - 2^2y iix = (17) 

x 2 + y 2 - 2uv - 2Hu 2 - V2(x 2 - y 2 )xu - 2\/2yx yu = e 

(18) 

Here, the overdot denotes a proper time derivative and 
H is: 

H = f + f + 2(x 4 - y 4 ) + 2 labX a x b + V2xv (19) 

We have considered the normalization condition for 
the 4- velocity, Eq. QB^ instead of the (complicated) 
equation for v. There, e = —1,0 for timelike and null 
geodesies, respectively. 

It can be advantageous to use u as affine parameter 
instead of the proper time. If we assume u > 0, from Eq. 



G3: 



1 d (hm) 
V2 du 



(20) 



Now, in the hyperplane y = 0, by inserting Eqs. (|20|) in 
(I16II . 1181) and recalling x = u (dx/du), a coupled system 



On the other hand, in the hyperplane x — Eqs. (|15(l 
and l|17l) simplify to 



u = u = constant = 7 



ii z H, y = 



(21) 
(22) 



taking the same form as in the type N reduction (see 
for example |13|). Analytical differences come from the 
quartic term of i|19|) which is absent in the type N. This 
difference can be physically illustrated when we consider 
the functions / = k(u) z 2 , with n(u) an arbitrary real 
function 4 , and \ a = A a (u) z in Eq. (Tg) . Then, the 
motion in the y direction is governed by: 

d 2 y/du 2 - 8y 3 + 2{ lab \ a {u)\ b (u) ~ n{u))y = (23) 

This equation describes a (parametric) non-linear oscil- 
lator in a potential V(y) = y 2 (j a bX a X b — k — 2y 2 ). In 
the type N reduction there is no cubic term in Eq. (|23|) 
and the motion reduces to that of an harmonic oscillator. 

Non-linear oscillators can have chaotic motion. 
Geodesic chaotic motion has been found in the re- 
lated case of non-homogeneous vacuum pp waves 
(/ ~ z n , n > 2), see |2l|. However, chaotic behavior is 
not always present in non-linear systems and therefore 
possible chaotic motion arising from Eq. (|23|l would 
need further investigation. 

• Case W, v = -2/(z + z) 

We take the u dependence out of W° respect to the 
function considered in subsection B, so now W° = z. 
Geodesies in this second case have the more complicated 
form: 

it- (1 - v/x 2 )u 2 - 2/x ux = (24) 

x + u 2 (H, x +{x - 2v/x)(\ - v/x 2 )) + 2 (l - 2v/x 2 ) ux 

- 2/x vii = (25) 

y + it 2 (H, v +y(l - v/x 2 )) + 2y/x ux = (26) 

x 2 + y 2 — 2uv — 2Hu 2 — 2(—2v/x + x) xu + 2y yu = e 

(27) 

Again, the normalization condition of the 4-velocity 
takes the place of the equation for v. The function H is 



H = H 



2x 2 



(28) 



and H° is given by l|14f) but with g(u) = 1. 



4 This choice is well motivated as it corresponds to homogeneous 
pp waves in the (vacuum) type N reduction. 
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The geodesic equations can again be rewritten as dif- 
ferential equations in terms of an affine variable u instead 
of the proper time. In particular, for the first equation: 

^ \ ( ^\ 1 V (29) 
du \ x 2 J x 2 

Then the remaining geodesies become equations for 
(x, y, v) in terms of u. 

The existence of a singularity at x = is clear from 
the geodesic equations, but not its character. Geodesies 
for a general vacuum type N or conformally flat pure 
radiation solutions have been analyzed in j22(, where 
it was concluded that x — is a physical singularity. 
Whether this is also the case for our type III solutions 
remains an open question. 



general solutions exist, as the two type III (classes of) 
EYM solutions presented here are also null and twistfree. 

It has been proved [24| that all spacetimes with 
vanishing curvature invariants of all orders: 1) are 
Kundt metrics of Petrov type III, N or O, and 2) are of 
Plebahski-Petrov (PP) type N or O (null radiation or 
vacuum). These results have been generalized to higher 
dimensions in (25j). Such spacetimes have been called 
VSI spacetimes. Our solutions clearly meet requirements 
1) and 2) and are therefore new examples of D = 4 
VSI spacetimes. As such, they might have interesting 
physical applications in supergravity and string theory, 
in the fashion of pp waves. 



III. CONCLUSIONS 

In [23 Tafel claimed that null twistfree solutions of 
the Yang-Mills equations were exhausted by the type N 
(class of) solutions given in [(|. A new type N solution 
was recently given in |l6| . We have shown that more 
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